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We construct a one-parametric family of the double-scaling limits in the hermitian 
matrix model 4>® for 2D quantum gravity. The known limit of Bresin, Marinari and 
Parisi belongs to this family. The family is represented by the Gurevich-Pitaevskii 
solution of the Korteveg-de Vries equation which describes the onset of nondissipative 
shock waves in media with small dispersion. Numerical simulation of the universal 
Gurevich-Pitaevskii solution is made. 

l.In this paper we construct a one-parametric family of the double-scaling limits 
in the hermitian matrix model 4>® for two-dimensional quantum gravity. This family 
is described by a common solution of the Korteweg-de Vries (KdV) equation 

T VVx T Vxxx 0 (1) 

and the following ordinary differential equation (ODE) 

Vxxxx + 5vVxx/3 + 5{vx)‘^/Q + 5{x -tv + n^)/18 = 0. (2) 
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The leading term of the asymtotics of this solution as |x |—^-cxd is the solution of the 
cubic equation 

X -tf + f = 0. (3) 

In the particular case f = 0 the solution of (@),® is reduced to the known solution 
of 0. On the other hand, the common solution of (|l|) and (H) exactly coincides 0 
with the well known Gurevich-Pitaevskii (G-P) special solution of the KdV equation 
0, 1^ which universally describes the onset of oscillating dissipationless shock waves 
in media with small dispersion. 

2.An essential development in the theory of two-dimentional quantum gravity 
based on the study of nxn hermitian matrix models 

Zn = J dHexp{-(3TrU{H)), (4) 

where if-nxn-hermitian matrix, U{z) = was made in the series of works 

The progress was achieved, considering the double-scaling limits in models 
(0), (h—0): n/l3={l + 5xh‘^)A. The following remarkable circumstance 

was then observed 0,0: 

It turns out 0,0 that the calculation of the second derivative of the limit of the 
non-regular part of In ( its regular part is irrelevant) is reduced to the finding of 
the limit of solutions of nonlinear difference equations 

(d Qi^Rni Rn+l) Rn—1) •••; Rn+U Rn-l), (5) 

which the right hand sides Q uniquely defined by the potentials U{z). It was es¬ 
tablished that for any U{z) there exist a natural m > 1 and constants A, B and 5, 
such that the asymptotics of solutions Rn of (0 is described by the formula (p is a 
constant): 


Rn = p{l + h^/^v{x)+ 

where functions v(x) satisfy of the first Painleve ODE ( for m = 2) and its higher 
analogues ( for m > 2). Before passing to a detailed discussion of a particular case 
of the potential 

U{z) = z^ 12 +giz^^ + g 2 Z^, (6) 
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which we consider in this article, we would like to give the following useful general 
statement which should be taken into account when calculating double-scaling limits: 
For any integer k the functions v{n + k) are expanded into Tailor series 

CXD 

v{n + k) = v{x + hk) = v{x) + '^{hkyv^^\x)/{l\), (7) 

i=i 


in which the constant c does not depend on the set of gj of (^), the natural m, and 
the constants A,B,j3,S,p. 

Remark. To make the results of the present paper compatible with Q and , we 
use notations different from []^ . This is the reason why we fix the following particular 
value of the constant c: c = \/ 6 . 

3.The equation for the potential @ has the form ||TI| : 


n/(3 = Rn{l + 4:Ji{n)+ 6 J 2 ), (8) 

Ji{r) Rn+l T Rn T Rn—ly ' 72 ('^) T Rn+ 2 Rn-\-l Rn+lRn —1 T Rn—lRn— 2 - 

In the case of the general position , one has the case m = 2. Assuming that the 
constants 6 , p,A,gi and g 2 do not depend of h, the substitution of 

n/j3 = A{1 + 6 h‘^x), Rn = p(l -l- h‘^v{x) + ...)/3 (9) 


into @ and equating the coefficients at different powers of h to the zero gives the 


following sequence of relations [11 


+ 12p5'i -h 1) = A, (10) 

■ p\2AQpg2 + 2Ag^) = -2A, ( 11 ) 

■ -{91 + 10pp2)^a: = {lhpg2 gi){2v^^ + n^). (12) 

It follows from the first two relations that p is a solution of the quadratic equation 


180p2P^ + 245 (ip 1 = 0. 


(13) 


Except for the degenerate case 


92 = -Pi/(15p) = 4p2/5, 


(14) 
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(i.e. the case of the multiplicity of roots in ([T3|)), the function v{x) turns out to be a 
solution of the hrst Painleve equation. 

In the degenerated case, one needs another limiting transition. Assuming that (|^) 
is true, Bresin, Marinari and Parisi considered in |l| the following limiting transition: 


Rn = -(1 + h‘^v{x))/12{gi),n/P = -(1 - h^x)/36{gi), 


where v is solution of ODE 


(16) 


Vxxxx + 5vVx^/3 + 5{v^Y/() + 5(x + ■u^)/18 ^ 0, 


V(X} 


|x| 


oo. 


(16) 

(17) 


The numerical solution of the boundary value problem (p!6D, (p!7D was obtained in 
( it was based on a difference scheme in which the ODE ([T6| ) was replaced by the 
starting discrete equation (^).) That numerical calculation has shown the uniqueness 
of the solution of the boundary value problem ©• 

The work |]^ had a wide response. The limiting transition considered in that paper 
was studied later from different points of view in a series of articles (see e.g. example 
jT^,|T^, ). An unexpected connection of that transition with the an old problem of 
the onset of dissipationless shock waves was discovered in 1994 in p| . It was observed 
that the solution of the boundary value problem (pTl) coincided with the G-P special 


solution of the equation (|^ at f = 0. 

However, in spite of the signihcance of |^, the given analysis of the limiting 
transition in the degenerated case was not satisfactory from the view point of a general 
requirement for investigations of degenerate cases of that kind. According to the same 
requerment going back to H.Puancare, ”... the investigation of degenerate systems 
should not be restricted by the study of the picture in the point of degeneration, 
but should include the description of the reorganizations which take place when the 
parameter passes through the degenerated value” . 

It will be shown in the next section if we try to satisfy that requirement and 
take into account the additional statement (0). we are led to the replacement of the 
boundary value problem ([l^ ) , (|T7| ) by the one-parametric family of boundary value 
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problems (||),(|^). The particular case of the double-scaling limit in the degenerate 
case ( 0 ) exactly corresponds to t = 0. 

4. Assume that the critical difference gi — has the order 0{h^)\ 

gi = a + bh^ + ...,g 2 = 4a^/5 -|- rbh^ + .... (18) 

(Here p,a,b and r are constants independent on h.) It follows from the square equation 
([T^) follows that in this case 

p = -l/(12a) + (19) 

and the relation (p!2D implies that the parameter b should also be small: 

5 = 7 / 1 " + ... (20) 

(^,7 are constants independent on K). Assuming as before that ®. (Q) are true, 
substituting the relations (0), (|^), ([l8|) - (^0D into (||), and then equating to zero the 
coefficiets at h\ we autoumatically (upon excluding of irrelevant trivial cases) come 
to the conclusion that 


p = 4,q = 2,gi = a(l - th'^/S + ...),g 2 = 4a^(l - + ...)/5. (21) 


The parameter t is calibrated in accordance with notations of |^. 

Taking into account the relations (i, (H) and (^), we obtain after some changes 
of notation a one-parametric family of double-scaling limits 

n/{3 = —(1 — h^x + ...)/(36a), = — (1 + h‘^v{t, x) + ...)/(12a), (22) 


where n is a solution of the ODE (^.(The last statement is obtained if one substitutes 
,(p^) into the discrete string equation (^ and equates coefficient at to zero.) 

5. The solution v{x) of the ODE (pTf ) of satisfies the boundary conditions (|^. 
On the other hand, the ODE (^ is isomonodromic |^, and the complete set of its 
monodromic data is completely defined by the boundary condition ([I7D.( It has been 
calculated at / = 0 by Moore in ) This implies (see @]) that the solution v{t,x) 
of the ODE d^) is at the same time the solution of the KdV equation (|I]), and that the 
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leading term of its asymptotics as |x| —>• cxo is the solution of the cubic equation (H). 
According to 0, exactly this feature defines a special solution of the KdV equation 
which describes in a general way the onset of dissipationless shock waves in problems 
with small dispersion. 

Actually, the corresponding solution (|^) is the leading term of the asymptotics 
of v{t,x) as t —> —oo for any x as well. If t —> oo, this asymptotics is true for any 
X except for a limited (but expanding with the growth of t) area filled with high 
frequency oscillations which correspond to the process of dissipationless shock wave 
generation. 

For the description of the behavior of n(t, x) in this area it was suggested in to 
use the self-similar solutions of the Whithem equations which arise after averaging 
with respect to the period of a ’’knoidal” wave (one-phase periodic solution of (Q)) 
The self-similar substitution suggested by Gurevich and Pitaevskii (correlating to the 
self-similarity of the solution of (^)) reduces the Whithem equations to the ODE. 
which was numerically investigated in Later a corresponding solution of the ODE 
was found in an explicit form by Potemin [^. ( This solution of Potemin, as has 
been pointed in |Q] , may be easily found, using the results of and the fact that the 
ODE (j^) is valid for n(f, x). ) Potemin also has found the exact values of the so-called 
trailing s_ and leading s+ fronts, which define the boundary of the oscillation region 
(s = h 

s_ = -72,8+ = 7 iO/ 27 
numerically established in before. 

The next step to clarify the behaviour of n(f, x) was done in the recent paper 
in which the uniform behavior of the GP solution in the neighborhood of the trailing 
edge was studied (the answer is obtained in the terms of the separatrix solution of the 
second Painleve solution). Besides, the speculations of 0 show that the ’’average” 
description of [^] and the real leading term of G-P solution asymptotics are probably 
the same within the oscillation region up to the shift of phase of 7r/2. 6. However, 
those results concern only to the behavior of n(f, x) with great t. Taking into account 
the universal nature of the G-P special solution, the problem of numerical simulation 
of the G-P solution’s behavior with a finite |f| (solved in the last section of this article) 
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is of interest as well. To solve this problem, we use an elementary iterative scheme 
similar to one of []^ in which ODE (^) is replaced by a discrete equation equivalent 
to the discrete string equation (g). Namely, the following procedure is used ([—L, L] 
is a numerically investigated interval, e = L/N is the step of the uniform net): 

Vn+l{k) = Vn{k) + P[Vn{k), fc], 

P[v{k), k] = Rk[l + 4:giJi{k) + Qg 2 J 2 {k)] - (1 - e^/c\/6)/3, 

= Rk+l + Rk + Rk-l: J2{k) = (Ji(fc))^ + Rk+2Rk+l ~ Rk+lRk-l + Rk-lRk-2, 
gi = -il- teV3)/12, g 2 = {I - t€^)/180, Rk = {l + e^Vk), 
no(0) = 0,r;o(^) = —{ekV6) ^ — t{ekV6) ^ /3{k ^ 0) 

and in all approximations for \k\ > N — 1: Vn{k)=VQ{k). Results of the numerical 
analysis 1, 2 correspond to the results of 0 on the qualitative level. In particular, a 
good correspondence has been observed between theoretical results and the positions 
of the trailing s_ and the leading s+ fronts of the oscillation region. 

Authors are grateful to V.E. Adler for the help in the numerical simulations. This 
work has been supported by the Russian Foundation of Fundamental Researches 
(96 — 01 — 00382). The second author was maintained by Foundation for Leading 
Scientihc Schools of Russia ( 96-15-96241). 
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